Combinatorial Bethe ansatz 
and ultradiscrete Riemann theta function 
with rational characteristics □ 



Atsuo Kuniba and Reiho Sakamoto 

Abstract. The U q {sl<i) vertex model at q — with periodic boundary 
condition is an integrable cellular automaton in one-dimension. By the 
combinatorial Bethe ansatz, the initial value problem is solved for arbi- 
trary states in terms of an ultradiscrete analogue of the Riemann theta 
function with rational characteristics. 

1. Introduction 

Solvable vertex models [5] associated with a quantum affine algebra U q yield 
completely integrable one-dimensional cellular automata at q = [9] [8] [6] . The 
simplest case of U q (sl2) with periodic boundary condition is known as the periodic 
box-ball system [TTl |2"2"] . Here is an example of the time evolution pattern: 

t = 0: 111222221111222222221111111111111112222111211 
t = 1: 111111112222111111112222222221111111111222122 
t = 2: 222211111111222211111111111112222222221111211 
t = 3: 111122222222111122221111111111111111112222122 
t = 4: 222211111111222211112222222221111111111111211 
t = 5: 111122221111111122221111111112222222221111121 
t = 6: 222211112222111111112222111111111111112222212 
t = 7: 111122221111222222221111222221111111111111121 
t = 8: 111111112222111111112222111112222222221111112 
t = 9: 222211111111222211111111222211111111112222221 
The dynamics is described as a motion of balls (letter 2) hopping to the right along 
the periodical array of boxes. The fusion transfer matrices T%,T2, ■ ■ ■ in the q = 
vertex model give rise to a commuting family of deterministic time evolutions. The 
above example is the evolution under Tg. 

In [17] , the initial value problem of the periodic box-ball system is solved by an 
inverse scattering method. It synthesizes the combinatorial versions of the Bethe 
ansatz [3] at q — 1 [13] and q = [14]. The action-angle variables are constructed 
from the rigged configurations (q = 1) by imposing a certain equivalence relations 
specified by the string center equation (q = 0). The time evolution T*(p) of any state 
p is obtained by an explicit algorithm whose computational steps are independent 
of t. 

The Bethe ansatz approach [17] captures several characteristic features in the 
quasi-periodic solutions of soliton equations [H |TT] . It generalizes the connection 
between the original box-ball system on infinite lattice [20] and soliton solutions 
via the ultradiscretization [21] to a periodic setting. For instance, the nonlinear 
dynamics of the periodic box-ball system becomes a straight motion of the Bethe 
roots (angle variable) which live in an ultradiscrete analogue (14. ip of the Jacobi 
variety. 
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In this paper we express the algorithmic solution of the initial value problem 
[T7] by an explicit formula involving the ultradiscretization of the Riemann theta 
function with rational characteristics (z G R 9 , a G (Q/Z) 9 ): 

6 a (z) = lim elogtf a (z) 

(1 1) 

v ' ; = - min{*(n + a)f2(n + a)/2 + t (n + a)z}. 

nez» 

Here $ a (z) is the Riemann theta function [19 with a pure imaginary period matrix: 

(1.2) * a (.) = £ ex P (- t(n + a) " (n + f 2 + t(n + a)z ). 

nGZs 

Here Q is a symmetric positive definite integer matrix (|5.ip . The function a (z) 
enjoys the quasi-periodicity analogous to the Riemann theta function: 

(1.3) 6 a (z + v) = 6 a (z) + t vrr 1 (z + v/2) forany veffl 9 . 

Our main Thcorcm l5.1l covers all the states in the periodic box-ball system, which 
generalizes the earlier result [15] on those states in which the amplitudes of solitons 
were assumed to be distinct. The theorem is derived from the recently obtained 
piecewise linear formula [TB] for the Kerov-Kirillov-Reshetikhin (KKR) bijection 
[13] from rigged configurations to highest paths. It involves the ultradiscretization 
of the tau function for the KP hierarchy [10 . A state p of the periodic box-ball 
system is effectively treated as the infinite system p®p®p®- ■ ■ , and the ultradiscrete 
tau function for such states turns out to be expressible in terms of O a (z) as in (|5.5[) . 

In section [5] we define the periodic box-ball system. In section [3] various data 
in the combinatorial Bethe ansatz at q = 1 [13] and q = [14] are introduced. In 
section [4] we recall the solution of the initial value problem by the inverse scattering 
transform. In section [5] we present our main Theorem 15. II in this paper. 

2. Periodic box-ball system 

Let us recall the periodic box-ball system without getting much into the crystal 
base theory [T2]. For a positive integer let B\ = {(21,2:2) G (Z> ) 2 x\ + 22 = 1} 
and set ui = (I, 0) E B\. The element (21, X2) will also be denoted by the sequence 

1 ... 1 2 ... 2, which is the semistandard Young tableau (without a frame) containing 
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the letter i xi times. In this notation ui = 1 . . . 1, and the two elements (1,0) and 
(0, 1) in B\ are simply written as 1 and 2, respectively. In the following, the symbol 
£g) meaning the tensor product of crystals can just be understood as a product of 
sets. Define the map R : Bi ® B\ — » B\ (g) B\ by 

'l<g(Z,0) if (x!,x 2 ) = (Z,0) 

1,22 — 1) otherwise, 



(21,22) ®1h i 

I 2 (gi (21 

(21, 22) g) 2 



J 2® (0,0 if (xi, £2) = (0,0 

1 1 (g) (21 — 1, 22 + 1) otherwise. 

i? is a bijection and called the combinatorial R. We write the relation R(u (g> 6) = 
&' g) «' simply as w ® 6 ~ b' g> w', and similarly for any consequent relation of the 
form a<g>u(g&<g>c~ag)6'g)u'cg)c. 
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A state of the periodic box-ball system is an array of 1 and 2, which is regarded 
as an element 61 eg) • • • Cg> bh G Bf L with L being the system size. Let the number 
of 2 G B\ appearing in b\ C§ • • • CS> &l be N. We assume L > 2N. The other case 
can be reduced to this case by the exchange 1 <-> 2 as explained in §3.3 of [17]. Let 
V be the set of such states. Then the time evolution Tj : V — > V is defined by 

(2.1) m®p~p*®vu vi®p~Ti(p)®vi. 

In the first relation, one applies the combinatorial R L times to carry ui through 
p G V to the right. It determines v\ G B\ and p* G V uniquely, (p* will play no role 
in the sequel.) Then the second relation using the so obtained vi specifies Ti(p), 
where the appearance of the same vi in the right hand side is a non-trivial claim 
( [17] (2.10)). vi is dependent on p as opposed to ui. 

The combinatorial R is the identity map on B\ eg) Bi, therefore T\ is just the 
cyclic shift Xi(6i<g>- ■■®b L ) = b L ®bi®-- -®b L -i- The commutativity T{F k = T k Ti 
holds for any k,l ([II] Theorem 2.2). 

Example 2.1. We depict the relation R(u <E> b) = b' eg) u' by a vertex diagram: 

b 



Then for p = 1121221 G Bf 7 (symbol ® omitted), one has v 2 = (1, 1) = 12 € B 2 
from (|2.ip . The time evolution T2(1121221) is obtained by the composition of the 
vertex diagrams 

112 12 2 1 
12+11+11+12+11+12+22+12 
2 112 112 

leading to T 2 (1121221) = 2112112. 

3. Combinatorial Bethe ansatz data 

Let 'P_|_ be the subset of V consisting of the s^-highest states. Namely, V+ = {p£ 
V I eip = 0}, where e\ is a Kashiwara operator [12]. In practice, p = b\® ■ ■ ■ ®b^ 
is highest if and only if b\ <8> • • • ® h% has nonnegative weights for any 1 < k < L. 

According to the combinatorial Bethe ansatz at q = 1 [13], there is one to 
one correspondence between V+ and the set of the objects called the rigged con- 
figurations. We let <fi(p+) = (/i, J) denote the Kerov-Kirillov-Rcshctikhin (KKR) 
bijection. Bere p+ is highest state and the right hand side is a rigged configuration. 
The data fi, called a configurations, is just a Young diagram for a partition of iV", 
for which N < L/2 is assumed. 

In what follows, we let I — {i\ < i 2 < ■ ■ ■ < i g } be the set of distinct lengths of 
rows of /z, and the multiplicity of the length i(G I) rows by m^(> 1). In terms of 
the non-increasing array of parts we have 

(3.1) / i=(zr 9 ...c i )- 

Obviously, a configuration (i is equivalent with the data I and , . . . , rrii g . 

To each row of ^, there is assigned an nonnegative integer called a rigging, which 
is represented by the symbol J in {[i, J). See Example 15. 21 Let J^i, . . . , Ji, mi be the 
ones assigned to the rrii rows of length i in fi from the bottom to the top. The data 
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[fx, J) with J — {(J;. Q ) | i G 1, 1 < a < to^} will be called a rigged configuration if 
< Jj,i < ■ ■ • < Ji.rm < Pi holds for each i e I. Here pi is defined by 

(3.2) pi = L - 2 min(z, j)m 3 -, 

which is nonnegative if /i is a configuration. 
We introduce the index set 

(3.3) I = {(«, a) | i e 1, 1 < a < mi}, 7 = |I| = m tl H h m ig 

labeling the rows of fi. The following vectors and matrix whose components are 
labeled with I will be used: 

(3.4) h J = (mm(i 1 j)) {l a)ei eZ-< (j £ I), 

(3.5) 3=(J iia + a-l) (i>a)ei eIP', 

(3.6) A = (Aa,jp)(i ia )^)^ > AiajP = 5i,j8a,l3(Pi + m i) + 2 j) ~ ^J- 

We also introduce the vectors and matrices whose indices range over the previously 
introduced set I = {ii, . . . , i g } as 



(3.7) 




(min(i, j)) ie j e Z 9 , 


(3.8) 


P = 










(3.9) 


J = 


+ • • • + J^mjiel S Z 9 , 


(3.10) 


F = 


(Fi,j)i,jei> Fi.j — ^ ] Aj a jp = S{j 






0=1 


(3.11) 


M = 


diag(m l ) ie i, 



The matrices A, F were introduced in the study of q = Bcthc equation [T3]. In 
particular, ^4 is the coefficient matrix of the linearized Bethe equation (string center 
equation), which is known to be positive definite. The result of the sum in (| 3 . 1 [) 
is independent of a. 

4. Action-angle variable and inverse scattering transform 

Any state p 6 V can be made highest by applying the cyclic shift appropriately. 
Namely, one can express p as p — Tf(p + ) with some integer d and highest state 
p + G V+- Although the choice of d and p+ is not unique, the configuration of 
4>(p+) is unique. The set of states is decomposed according to the configuration of 
p + as V — U Jli 'P(/i). We call fi the action variable. It is a conserved quantity f [IT] 
Lemma C.3, Corollary 3.5), therefore each subset V(fJ.) is invariant under any time 
evolution Tj. The action variable fi is equivalent with the data I = . . . , i g } and 
{rrii | i G 1} via (|3 . 1[) . In the periodic box-ball system, I is the list of amplitudes of 
solitons and m; is the number of solitons with amplitude i. 

Given an action variable fj,, the set of angle variable is given by 

(4.1) J{ji) = (l mii x ■ ■ • x X mig )/r, Y = AP, 

where I n = (Z n — A n )/S„ is the n-dimensional lattice without the "diagonal 
points" A„ = {(zi, . . . , z n ) £ Z™ | z a = zp for some 1 < a ^ (5 < n} divided by 
the action of the symmetric group & n . Plainly, it is the set of array of n distinct 
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integers whose ordering does not matter. In the sequel, elements in I mil X ■ • ■ xl mij 
and its image in will be denoted by the same symbol. The set J{[i) is finite 

and its cardinality is given by [HI [T7] 

(4 „ lwl . (detF) n±(-:--)^n(-;r 

This is a positive integer if n is a configuration. (When = L/2, we have pi = 0. 
In this case, the factor — ( Pi s +mi s _1 ) is to be understood as -^—.) The formula 
()4.2p tells the number of states containing im solitons with amplitude i. The 



completeness of soliton states is known in the sense that ^1^(^)1 = (at) f° r 
< N < L/2, where the sum extends over those such that J2i^ m i = N ([14j. 
Theorems 3.5 and 4.9). The states in section Q] have the length L — 45 and fi 
depicted in Example l5.21 For this fi and L, the formula (|4.2[) tells \J(n)\ = 316350. 

Let us specify the angle variable for a given state p 6 (/•*)• Suppose p is 
expressed as p — Tf(p + ) in terms of a highest state p+ and the cyclic shift T\. Let 
the image of the KKR bijection be <j>(p+) — (fJ>, J). From the configuration /i we 
know I = {ii < %2 < ■ ■ ■ < i g } and , . . . , rrii g by (13 - 1 j> - Then the angle variable 
is found by applying the definitions (j3.4|) . (|3.5p . (|3.6p to J = (Jj, a )(j Q ) e j as 

$ : 7>(/i) — » ZxP + ^ 

(4.3) 

p ► (d,p+) ► (J + dhi)/r. 

Due to the shift +a— 1 in (|3 . 5|) and the rigging condition < J^i < • • ■ < Ji, mi < Pi, 
it follows that the components of J hence J + dhi within each block labeled by i e I 
are distinct. In (|4.3[) we regard the J + dhi as a representative of an element in 
j7(/z). The non- uniqueness of d, p+ is cancelled by modT making $ well-defined 
( [T7] Proposition 3.7). 

Theorem 4.1 ( |17j Theorem 3.12). 3> is a bijection and the following diagram is 
commutative. 

VfjJL) — J(/i) 
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— 5— J( M ) 

//ere i/ie izme evolution on J(p) is defined by Tj(J) = J + h;. 



The composition o TJ o $ achieves the solution of the initial value problem 
by the inverse scattering transform. Namely, $ ±1 is the direct/inverse scattering 
map which transforms the dynamics T; of the periodic box-ball system into the 
straight motion on J'(fi) with the velocity h;. In this sense, serves as an 

ultradiscrete analogue of the Jacobi variety in quasi-periodic solutions to soliton 
equations [H I5"l [Tl"] . 

By using Theorem 14. 1[ a closed formula for the fundamental period, i.e., the 
smallest positive integer such that Tj^(p) — p for any I and p has been obtained in 
Theorem 4.9 of [T7], which supplements a recursive description of the I — oo case 
|22j . The states in section [T] have the fundamental period 3515. See also [TB] for a 
generalization to A n based on the Bethe eigenvalue at q = 0. 
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We note that the time evolution of the angle variable J (|3.5|) leads to the rule 
(4.4) T,(J)=J + Afh, 

for the "bundled" angle variable J (|3.9[) . where M is defined by p. lip . 
5. Explicit formula for the solution of the initial value problem 
Fix the action variable fi (|3.1j) and set 

(5.1) tt = MF = (f2i,j)i,jel) Ojj = <\jl>,n>, + 2mm(iJ)m l m v 

where F is defined in (13 . 1 1 [) . Q is a positive definite symmetric integer matrix which 
will play the role of the period matrix. Note that 

(5.2) f!hi = LMhi. 

In addition to the ultradiscrete Riemann theta function a (z) in we need 

x(s;I). It is a function of s = (s,) ie i £ Z 9 as well as I = (ij lQ + a — 1)^ a , £ | G 
X mii x • • • x I,„ is that takes rational values: 

(5.3) x(s; I) = x(s'; I) if s ee s' mod MZ 9 , 

(5.4) x( S ;i) = E^: E (^-^-f) 

if < Si < rrii, < • • • < li.mj for all i e I. 

Here the second sum in (|5.4[) extends over the pairs (a, /?). In particular, x(0; I) = 0. 
The condition J^i < • •• < ii, mi is needed to cope with the & mi symmetry in T mi . 

Given I G Z 9 and I = [Ji. a + a — 1)^ Q j £ j G X m<i x • ■ ■ x , we define the 
ultradiscrete tau function by 

(5.5) r r (k)= max {9 M - ls fl + Mtrh^ - kin - + X (s; I)} 

for r G {0, 1}, fc G Z. Here h,, p and M are those in ([3J]) . (|3T8]) and (fBTTT]) . 

Let us specify I and I in (|5.5[) from a given state p G V^i). As for I, it is taken 
to be the angle variable of p e 7 : '(^) according to (|4.3p . and I is determined from it 
as follows: 

(5-6) I = Hp) = (Ii, a +a- a)6i G J{n), 

(5.7) I = (/ <tl + • • • + h^iei G Z»/FZ». 

This is the same relation as (|3.5j) and (|3.9| . Note that the elements I and I are 
defined mod T = Al? and mod FZ 9 , respectively, which is consistent with the 
relation (|3.10p between A and F. Now we state the main result of this paper. 

Theorem 5.1. The state p G V(n) corresponding to the action variable /i and the 
angle variable I is expressed as p = (1 — x\, x\) <S> • ■ ■ <8> (1 — xl, xl), where Xk = 0, 1 
is given by 

(5.8) x k = r (fc) - r (fc - 1) - n(fc) +n{k- 1). 

The proof is similar to section 3 in 15J. Combined with Theorem 14. 1\ Theorem 
15.11 leads to an explicit formula for the solution of the initial value problem. The 
state (Y[ l Tf')(p) is obtained just by replacing the angle variable as 

(5.9) Ii-fi + ^qh, 

i 

6 



in (|5.5p . Under the change (|5.9p . the function x( s ', I) l|5.4|) is invariant since it only 
depends on the difference Ii t p — Ii_ a . On the other hand, the variable I entering 
(15. 5|) changes as 

(5.10) Ih-vI + MCj^qh,) 

according to the rule ijO]) . Especially, T/^I) = I + LMhi = I + fihi by ([572]) . 
Thus the tau function only exhibits the change due to the quasi-periodicity (|1.3|) 
leaving (|5.8|) invariant. This is consistent with the fact that the cyclic shift T± is 
trivial and (15. 8p is a period L function of fc. We remark that our tau function (|5.5j) 
is always integer- valued. 

From (|3.7p we see h; = hoo if I > maxl(= i g ). Therefore the effect f|5 . 10[) of 
the time evolution Ti is I i— » I + Mhoo , which is equivalent to setting r i— > r + 1 
in ()5.5p . In other words, r serves as the time variable for such 7], i.e., T*(p) = 
(1 - xi )t , £i,t) ® • • • ® (1 - x L , t , x L ,t) with 

(5.11) x ktt = n(k) - T t {k - 1) - 7-t+i(fc) + r t+1 (fc - 1) 
holds for I > max I. 

Example 5.2. Let us consider the t = 0,1 states in section [TJ 

p = 111222221111222222221111111111111112222111211 
^■ 12 ' ) T 9 (p) = 111111112222111111112222222221111111111222122, 

which are elements in Bf 45 . We have p = T? 9 (p + ) with a highest state 

p+ = 111111222211121111122222111122222222111111111 
which corresponds to the rigged configuration: 





















6 










2 




10 









Thus I = {1,4,9}, M = diag(mi, r?i4, mg) = diag(l,2,l). Ii iCt in (15. 6p reads 
(Ii,i,/ 4 ,i,/4,2,/9,i) = (39,31,35,29). The data p (|3Tgj) . I ((Ell) and fi ((53) are 
given by 

/pA /37\ /39\ /39 4 2 

p = p 4 = 19 , I = 66 , = 4 70 16 

Vp 9 / V 9 / V 29 / V 2 16 27 / 

The argument of the tau function (|5.5p reads 

/r-fc + 41/2\ 

(5.13) I + M(rhoo - jfehi - -) = 8r - 2fc + 47 =: z nfe . 

2 \9r-/c + 49/2y 

The function x(s;I) in (|5.4p is non-vanishing at s = (0,1,0) taking the value 
^-(-^4,2 — ^4,i — ^ ) = Thus the tau function is given by 

T r (k) = max{eo,0.o(Zr,fc)> B Q l n( z r,fc) - ^j-}- 

' 2 4 



The time evolution pattern in section [T] is reproduced according to (|5.11|) . Here is 
the plot of T r (k) for < r < 2 and < k < 45. 



T r (k) 
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The boundary values are (r (0), Ti(0), t 2 (0)) = (15, 28, 41) and (r (45), ti(45), t 2 (45)) = 
(13,8,3). In the plot of r t (fc) or r t+ i(fc), one observes a change of the gradient at 
those k corresponding to the front or tail of solitons in Tg(p) (|5.12|) . 

The tau function (|5.5p simplifies considerably if there are no solitons with the 
same amplitude, namely, Vrrii = 1. We have x( s ;I) = 0, I = I and D = F = 
{SijPi + 2 min(i, j))ijei- Let us write the ultradiscrete Riemann theta function 
with characteristics as O(z) = &o(z). Then r r (k) = 0(1 — f — &hi + rh^. 
Consequently (|5.8[) reduces to 

Xk = e(i- 1 - fch!) - e(i - 1 - (fe-i)hi) 
- e(i- 1 - fchx + hoo) + e(i- 1 - (fc-i)hi + hoo), 

which was obtained in |15j . This has the same structure as the ultradiscretization 
of the solution of the periodic Toda lattice by Date and Tanaka [4] . 
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